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Many modern theories predict that the fundamental constants depend on time, position or 
the local density of matter. Here we develop a spectroscopic method for pulsed beams of cold 
molecules, and use it to measure the frequencies of microwave transitions in CH with 
accuracy down to 3 Hz. By comparing these frequencies with those measured from sources of 
CH in the Milky Way, we test the hypothesis that fundamental constants may differ between 
the high- and low-density environments of the Earth and the interstellar medium. For the fine 
structure constant we find Aa/a = (0.3 ±1.1) x10 -7 , the strongest limit to date on such a 
variation of a. For the electron-to-proton mass ratio we find = ( - 0.7 ± 2.2) x 10 ~ 7 . 
We suggest how dedicated astrophysical measurements can improve these constraints 
further and can also constrain temporal variation of the constants. 
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n the Standard Model of particle physics, the fundamental 
constants are fixed parameters, but in many extensions of the 
standard model they may change 1,2 . In higher- dimensional 
theories for example, which aim at a unified description of gravity 
and the other forces, the constants depend on the size of the 
compactified extra dimensions. Their size may be changing in a 
conceptually similar way to the observed expansion of the 
Universe, causing the constants to change with position or time. 
Some theories of dark energy also predict a variation of the 
constants. These theories hypothesize that a scalar field is 
responsible for the observed cosmic acceleration. In chameleon 
models, this field acquires an effective mass that depends strongly 
on the local matter density 3,4 . In our high-density environment, 
the field has such a short range that it goes undetected in tests of 
the equivalence principle and searches for fifth forces. In these 
models, particle masses and coupling constants may also depend 
on matter density 5 . Given the potential role of these theories in 
answering outstanding questions in physics, it is important to test 
the prediction of varying constants. 

Atomic and molecular transition frequencies provide a 
natural way to look for variation of the fine structure constant, 
a, and the electron-proton mass ratio, fi. Their present-day time 
derivatives are strongly constrained by comparisons of atomic 
clocks 6,7 , while observations of high red-shift objects test 
changes on cosmological scales. Some analyses of optical 
absorption spectra at high red-shifts suggest a variation of a 
across the Universe 8 , but others find no variation at the level of a 
few parts-per-million 9,10 . High red-shift measurements using 
microwave and mm-wave molecular transitions, which can be 
far more sensitive to changes in the constants than optical 
transitions, give null results for Aa/a at levels between 8 and 0.8 
parts-per-million 11-14 . For Afi/fi, a recent study of methanol 
absorption lines at high red-shift found no variation at the 0.1 
parts-per-million level 15 . 

Measurements within our galaxy test whether constants 
depend on the local density. For changes in a with matter 
density the lcr limit is |Aa/a|<2 x 10 ~ 7 (ref. 16). For changes 
in fi a study of methanol maser lines sets a lcr upper limit, 
|A^|<2.9 x 10~ 8 (ref. 17). However, Levshakov et a/. 18,19 , 
comparing terrestrial and astrophysical microwave transitions in 
ammonia and other molecules, find an eight-standard-deviation 
systematic difference. This suggests a fractional change in ja of 
2.6 x 10 _ 8 when going from the earth to the interstellar medium, 
tentatively supporting the chameleon hypothesis. This intriguing 
situation calls out for new, high-precision measurements of the 
most sensitive molecular transition frequencies and for higher- 
quality astronomical observations. 

Here, we measure the frequencies of the 3.3 and 0.7 GHz 
A-doublet transitions of CH, which are very sensitive to changes 
in a and /u (refs 20,21). With these measurements and existing 
astronomical observations of CH in the interstellar medium, we 
constrain variation of the constants in our Galaxy, and open the 
way for future radio-astronomy measurements of higher accuracy 
to test more stringently whether the constants vary. Our 
spectroscopic method has Hz-level accuracy, and takes 
advantage of the short pulse duration and small velocity spread 
of cold molecular beams to control the interaction between 
molecules and an applied microwave field. This method can also 
be applied to other molecules of interest in chemical physics, 
metrology, particle physics and cosmology 22-26 . 



Results 

The CH molecule. Figure 1 shows the relevant energy levels in 
the ground electronic and vibrational state of CH, X 2 II(v = 0). 
These eigenstates of the Hamiltonian are also eigenstates of the 
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Figure 1 | Relevant energy levels in CH. Approximate frequencies are 
given in MHz and the sensitivity coefficients for the two A-doublet 
transitions are shown. 



parity operator, P, and of the squared total angular momentum 
operator, F 2 . In addition, they are approximate eigenstates of N 2 
and J 2 , where J = N + S and N = L + R. Here, S and L are the 
electronic spin and orbital angular momentum operators, and R 
is the rotational angular momentum operator. The spin-orbit 
interaction splits the N = 1 level into a pair of levels labelled by 
J — 111 and J — 312. Each is split by a Coriolis interaction into a 
A-doublet whose components are the parity eigenstates, 
\p= ±1> = |A = +l> ±(-l) J ~ s \A = -1>, where A is the 
quantum number for the projection of L onto the internuclear 
axis. It is these A-doublet frequencies that we measure. The 
interaction with the hydrogen nuclear spin, 1— 1/2, produces a 
hyperfine structure in each of the A-doublet components. The 
levels are labelled as (JP, F). 

These two A-doublet transitions have been observed towards 
numerous clouds in the interstellar medium 27-36 , and the 
3.3-GHz transition has also been observed in other galaxies 37 ' 38 . 
Figure 1 gives the sensitivity coefficients of the transitions to 
changes in the constants, K a and (refs 20,21), defined such that 
the rest frequency emitted by molecules in the astronomical 
source is co^ (1 + K a Aa/a + K^Afi/ fi), where co^b is the 
laboratory frequency. As these coefficients are large and 
different, a very sensitive test is possible when the transitions 
are observed together in the same gas. Until now, this test was 
hindered by a lack of precise laboratory frequencies. 



Frequency measurements. Using the apparatus shown in Fig. 2, 
we measure the laboratory frequencies by the Ramsey method of 
separated oscillating fields 39 ' 40 . The source produces short, cold 
pulses of CH molecules that pass through a state selector, a 
transmission line resonator and then a detector. The state selector 
preferentially populates one of the two parity eigenstates of the 
A-doublet being measured. Inside the transmission line resonator 
the molecules interact with a standing-wave microwave field. A 
first microwave pulse (a nil pulse), of duration t = 15jis and 
angular frequency cl>, is applied when the molecules are at 
antinode m : of the field. This prepares an equal superposition of 
the two A-doublet components, which then evolves freely for a 
time T at the angular frequency co 0 > developing a phase difference 
of ST relative to the microwave oscillator, where S — (co — a> 0 ). 
We choose the free evolution time to be T— ml/(2v 0 ) — t, where 
v 0 is the mean speed, X is the wavelength in the transmission line 
and m is an integer, so that the molecules are now at antinode 
mi + m. A second microwave pulse of equal duration and 
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Figure 2 | Experimental setup for measuring the A-doublet frequencies. A supersonic CH beam is produced by photo-dissociation of CHBr 3 . The initial 
state is selected either by optical pumping (for J = 1/2) or by driving the (1/2 ~ ) - (3/2 + ) transition with a mm-wave beam (for J = 3/2). The 
molecules travel through a microwave transmission line where the A-doublet transitions are driven, and are finally detected by laser induced fluorescence. 
Not to scale. 
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Figure 3 | Ramsey interference data. Dots: measured population in (1/2 ~, 1) as a function of the microwave frequency using the method of 
separated oscillating fields. Line: fit to the data using the model discussed in the text, (a) Wide frequency scan with a free evolution time of T=443|is. 
(b) Narrower frequency scans for three different free evolution times. Green: 458 (is. Blue: 380 (is. Red: 302 (is. 



amplitude then completes the population transfer with a 
probability 



x [Xcos(f) cos(^) - 2cksin(f) sin(^)] 2 . 

Here, X = \/ n 1 + 4(5 2 t 2 and ft represents a possible phase shift 
of the microwave field (modulo 2n) because of the change in 
position of the molecules between the two pulses. The 
population in one of the parity states is measured, as a function 
of co, by time-resolved laser-induced fluorescence. From the 
time -of- flight profile we select molecules with arrival times in 
the range t 0 ± 5t, where t 0 is the most probable arrival time, and 
dt~0.02t 0 . This ensures that all molecules used in the analysis 
are sufficiently well localized near the antinodes when the pulses 
are applied. 

Figure 3a shows data measured over a wide frequency range 
around the (l/2 + ,l)-(l/2 _ ,l) resonance. The line is a fit to the 
model b + aP(S), where b and a are constants and P(S) is given 
by equation (1). We set /3 = 0, and t and T to the values used in 
the experiment, leaving only the offset b, the amplitude a (which 
is negative when m is odd) and the resonance angular frequency 
co 0 as fit parameters. The model fits the data well and we use it to 
determine a> 0 . Figure 3b shows narrower frequency scans for the 
same transition and for three different values of T chosen to 



make m = 4, 5 and 6. Here, as in all our data, we see no 
dependence of co 0 on m, so we average together measurements 
taken using different values of m. We also see no variation of co 0 
when we subdivide the data according to the arrival time of the 
molecules. 

Our model assumes that ft is zero, as for a perfect standing 
wave. As the standing wave is not perfect, /3 is not exactly zero 
and there will be a systematic frequency shift of 3f= p/(2nT). Let 
the amplitudes of the two counter-propagating waves be 
proportional to (1±A), so that A is the amplitude imbalance 
between them. The microwave electric field is 



E = AyJ cos 2 (kz) + A 2 sin 2 (£;z) cos(co£ - 0), (2) 

where A is an amplitude, k = 2nlk and (modulo ri) 

(j) = tan _1 (Atan(/cz)). (3) 

When A is small, (j) changes very slowly around the antinodes. 
Suppose the first pulse starts when the molecules are at Z\ — z 0 + 
dz 0 , where z 0 is the position of an antinode, and the second pulse 
starts when the molecules are at z 2 = Z\ + (1 + e)mXI2. Expanding 
the trigonometric functions for small displacements around 
the antinodes (that is, Sz 0 «l, e«l), gives df~ Ae(v 0 //l), 
independent of Sz 0 to first order. For the J—112 measurement, 
taking A = 0.1 and e = 0.05 (see Methods), we would expect a 
velocity- dependent frequency shift of 0.06 Hz(ms -1 ) -1 - 
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A second systematic frequency shift is expected because of the 
motion of the molecules during the two short n/2 pulses. 
Including the Doppler shift in the expression for the lineshape in 
the case of perfect n/2 pulses and expanding about the resonance 
frequency, we find that a travelling wave results in the systematic 
frequency shift Sf D = (1 — 4/7t)(v 0 M)(t/T). For the counter- 
propagating waves this is further reduced by the imbalance 
factor A. For the J—l/2 measurement, with A = 0.1, t=15jis 
and T— 450 jis, we expect a shift of — 0.01 Hz (m s ~ x ) ~ l . 

We control these and any other velocity- dependent shifts, by 
measuring all the frequencies for at least three different values of 
v 0 . We find that the measured frequencies shift linearly with v 0 , 
and that the gradient, d/7dv 0 , differs for different values of 
The largest gradients we measure are 0.07 ± 0.01 Hz (m s _ l ) ~ 1 
for the /=l/2 measurement, and 0.03 ± 0.01 Hz (ms~ l ) ~ 1 for 
the J — 312 measurement. After extrapolating to zero velocity the 
results for different values of m x all agree, so we take the weighted 
mean of the zero -velocity values as the final frequency. This 
procedure does not rely on a knowledge of A. 

To minimize systematic shifts due to magnetic fields, the 
interaction region is magnetically shielded. For a linearly 
polarized microwave field and a small magnetic field, there is 
a symmetry in the shifts and amplitudes of the Zeeman 
subcomponents, and therefore no systematic shift. Circular 
polarization components of the microwave field are strongly 
suppressed inside the transmission line. For the 7=1/2 level the 
spin and orbital angular momentum contributions to 
the magnetic moment very nearly cancel, so Zeeman shifts are 
small. For example, the shift of the M F = + 1 component of the 
(l/2 + ,l) level is approximately O.OlHznT -1 . For the 7=3/2 
level the Zeeman shifts are far greater, 15.1M F HznT -1 
for F=l and 9.07M F HznT -1 for F — 2. To measure the 
residual field in the interaction region, we applied single 780 -jis- 
long pulses to molecules moving at 570ms -1 , and measured 
the Zeeman splittings of the J — 3/2 levels as a function 
of applied magnetic field in all three directions. Reversal of the 
applied field measures the residual field in that direction, 
averaged over the interaction region. We also look for any 
broadening of the single-pulse lineshape at zero applied field. 
The residual fields are consistent with zero in all three 
directions, with upper limits of 3, 56 and 25 nT along x, y and 
z, respectively. To measure the size of magnetic -field related 
shifts in the Ramsey data, we apply small homogeneous fields in 
each direction. For the J—112 transitions, we observe no 
significant shift at the 1 Hz level, even for fields as large as 50 uT. 
For the 7=3/2, AF = 0 transitions, the field-induced frequency 
shifts are all smaller than 0.04 HznT - l . Multiplying this by the 
largest upper limit to the field gives a systematic uncertainty of 
2 Hz. For the 7=3/2, AF= 1 transitions, we could not rule out 
gradients as large as 0.1 HznT - 1 for x, 0.2 HznT - 1 for y and 
0.05 Hz nT - 1 for z. Multiplying these by the upper limits to the 
residual field gives a systematic uncertainty of 1 1 Hz for these 
transitions. 

We measured the Stark shift of the J—112 transition by 
applying electric fields up to 20 V cm - 1 between the plates of the 
transmission line. Measuring the change in this shift upon 
reversing the applied field sets an upper limit to uncontrolled 
Stark shifts of 0.1 Hz. We checked for systematic shifts depending 
on microwave power or probe laser detuning, but these were 
negligible. Frequency shifts due to collisions in the beam, 
spurious frequency sidebands, blackbody radiation 41 , the 
motional Stark shift and the second-order Doppler shift are all 
negligible at the current accuracy level of 10 . 

Table 1 shows our measured frequencies along with their 
uncertainties, which are obtained by adding the statistical 
and systematic uncertainties in quadrature. The systematic 



Table 1 | Measured A-doublet frequencies with 1<r 


uncertainties. 




TVatiQitinn 




(1/2 + ,1)-(1/2-,1) 


3,335,479,356 ± 3 


a/2 + ,0)-(1/2-,1) 


3,349,192,556 ±3 


(1/2 + ,1)-(1/2-,0) 


3,263,793,447 ± 3 


(3/2 + ,2) -(3/2-, 2) 


701,677,682 ±6 


(3/2 + ,1)-(3/2-,1) 


724,788,315 ±16 


(3/2 + ,1)-(3/2-,2) 


703,978,340 ± 21 


(3/2 + ,2)- (3/2 "/I) 


722,487,624 ±16 


Levels are labelled with the notation (f, F). 



uncertainty is negligible for all but the last two, where it is 
11 Hz. The three J—112 measurements are consistent with the 
previous best laboratory 42 and astrophysical 30 ' 31 measurements, 
and are 300 times more precise than these. We know of no 
previous laboratory measurements of the J— 312 frequencies. 
There are astrophysical measurements of the first two J— 312 
transitions in the table 36 , and our results are consistent with these 
and about 300 times more precise. 

Constraining changes of a and \i in the Milky Way. With these 
new measurements and existing astronomical data, we can set 
useful bounds on the possible variations of a and \i with matter 
density. We study five interstellar sources where the density is 
typically 10 19 times smaller than on Earth. It is necessary to 
compare two spectral lines from each source, so that the Doppler 
shift can be removed, and for high sensitivity to changes in the 
constants these two should have very different sensitivity coeffi- 
cients. In one of the sources the J—112 and J— 312 A-doublet 
transitions of CH are observed together, and in the other four the 
J—112 CH transition is observed together with the ground state 
A-doublet transition of OH. Astronomical spectra are plotted on 
a velocity scale that is determined from the Doppler shift of the 
received radiation. For each selected source the spectral lines are 
narrow and have similar shapes, and we assume that the mole- 
cules producing them have the same velocities. To obtain our 
bounds, we assume that there is either a change in a or ja, but not 
in both (this is a limitation of our analysis). Then, the fractional 
change in a is 



Aa 

a 



1 



Avn + c 



CO, 



(1) 

lab 



(2) 
J lab 



(4) 



with a similar equation for the change in ji. Here, Av 12 — 
v (l) _ v (2) i s fae difference between the measured velocities of 
lines 1 and 2 and is determined from the published spectra, K$ is 
the sensitivity coefficient for transition i, co^ h is the frequency 
measured in the laboratory and co^ ro is the frequency used to 
define the velocity scale in the astronomical measurement. This is 
usually the accepted value for the transition frequency at the time 
of the measurement, is specified along with the published spectra, 
and carries no uncertainty. The laboratory measurements pre- 
sented here allow the second term in square brackets to be 
determined to high accuracy, so the uncertainty in Aa/a is now 
entirely due to the uncertainty in Av 12 . 

Table 2 shows the results obtained from the five selected 
interstellar sources. We find the central velocities by fitting 
Gaussian profiles to the spectral lines. Where there are multiple 
velocity components, we fit multiple Gaussians. In each source 
there are multiple velocity components and/or multiple hyperfine 
components, so we obtain several different values for the velocity 
difference. Their standard deviation is typically several times the 
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Table 2 | Analysis of astronomical data. 


Source Transition 1 Transition 2 


V 


AV12 


Av 12 


M (io ') 


(10 ') Ref. 




(kms 1 ) 


(kms 1 ) 


(kms 1 ) 






G111.7 - 2.1(CasA) CH(3264, 3335, 3349) OH(1667) 


-1.4, 0 


-0.01 (0.09) 


-0.08 (0.11) 


1.5 (2.0) 


-3.1 (4.1) 30,32 


G265.1 +1.5CRCW36) CH(3264, 3335) OH(1612, 1665, 1667, 1721) 


6.8 


0.06 (0.19) 


0.04 (0.16) 


0.9 (3.1) 


1.9 (6.4) 34 


G174.3-13.4(Heiles2) CH(3264, 3335, 3349) OH(1612, 1665, 1667, 1721) 


5.8 


0.00 (0.19) 


-0.02 (0.19) 


0.6 (3.6) 


-1.2 (7.4) 32,43 


G6.0 + 36.7(L134N) CH(3264, 3335, 3349) OH(1665, 1667) 


2.5 


0.05 (0.13) 


-0.12 (0.13) 


2.3 (2.4) 


-4.8 (5.0) 32,43 


G49.5 - 0.4(W51) CH(702) CH(3264, 3335, 3349) 


65 


-0.85 (0.53) 


-0.48 (0.55) 


-1.8 (2.0) 


3.6 (4.1) 36 


The columns are: the direction of the source, the transitions used (labelled by the frequencies of the hyperfine components in MHz), the velocity component(s) used, the measured velocity difference, the 


velocity difference corrected for any differences between the nominal and laboratory frequencies, the derived values of Aa/a and A/i/fi, the references where the spectra are given, 1(7 uncertainties are 


given in parentheses. 













individual uncertainties of the fits, probably because the 
underlying velocity distributions are not Gaussian. For each 
source, we use the mean and standard deviation of the 
individual results to give a central value and uncertainty 
for Av 12 . The angular size of the observed part of the source 
(the beam size) is often different for the two species. This can 
result in systematic velocity differences. This is not significant 
for CasA or W51, where the sources are smaller or comparable 
in size to the smallest beam size used. Neither is it significant for 
Heiles cloud 2 or L134N, because the beam size used for the 
CH and OH observations were nearly the same, and maps of 
these clouds show that velocity variation is less than 0.2 km s ~ 1 
over a much wider angular scale than this beam size 30 . For 
the observations of RCW 36, which is an extended source, the 
beam size for OH was twice that for CH. However, maps of 
this cloud 44 show that the observed region is centred on an 
intense localized clump whose size is approximately the 
OH beam size, and that any systematic velocity differences 
due to different beam sizes will be smaller than our quoted 
uncertainty. 

Using Av 12 , the measured values of co^ h for CH and OH 45 and 
the quoted values of co^ ro , we calculate the quantity in square 
brackets in equation (4), which we label as Av ' 12 in Table 2. From 
these, and the calculated sensitivity coefficients for OH (ref. 20) 
and CH (ref. 21), we find values for Aa/a and Afi/fi. Taking the 
weighted mean of their values, we arrive at the results Aa/a = 
(0.3 ± 1.1) x 10" 7 and A i u/ i u = ( - 0.7 ± 2.2) x 10 " 7 , giving us 
lcr upper bounds of |Aa/a| < 1.4 x 10 ~ 7 and |A^/^| 
<2.9 x 10~ 7 . 

Discussion 

Our analysis demonstrates the high sensitivity of this method. 
With dedicated astronomical measurements, far higher preci- 
sion could be reached and the study could be extended to 
sources at high red- shift. The scope is particularly great for 
studying sources where the 3.3 and 0.7 GHz transitions are 
observed simultaneously. There is currently a shortage of such 
data, and the comparison provides excellent sensitivity 
to variation of the constants and is relatively immune to 
systematic shifts that can arise when analysing transitions in 
different species. The relatively low intensity of the CH lines 
requires the use of large telescopes, but there are several 
that are suitable. The 100-m telescope of the Max-Planck- 
Gesellschaft covers the entire frequency range between the two 
A-doublets, whereas the 100-m Green Bank Telescope covers 
the frequency range from 0.29 to 2.6 GHz and so is suitable for 
measuring both transitions at moderate red-shifts. For a higher 
collecting surface, the 300-m telescope at Arecibo can be used. 
It covers frequencies from 3.3 GHz down to 1.225 GHz 
and selected frequency bands between 690 and 312 MHz, so 
could be used to study both transitions over a wide range of 
red-shifts. In a few years, even higher sensitivity will be 



provided by a new 500 m telescope being built in the Guizhou 
province of China. To circumvent problems of rf interference, 
interferometric measurements are efficient. The Very Large 
Array provides such measurements with high sensitivity 
and currently covers all frequencies between 1 and 3.3 GHz. 
Looking further ahead, the Square Kilometre Array, to 
be constructed in South Africa and Australia, will provide 
unprecedented sensitivity. Being an interferometric facility 
located in areas of low interference, this telescope should 
be able to measure the two A-doublets of CH in a large 
number of interstellar clouds, both at low and high red- shift. 
With a reasonable precision of 0.1 km s _ 1 in the measurement 
of Av 12 , and a modest dataset, it will be possible to 
reach below the 10 ~ 8 level for both Aa/a and Afi/fi. The 
laboratory frequencies are now so accurate that their 
uncertainties are unlikely to be a limiting factor in any future 
measurements. 

Methods 

Frequency measurements. The setup for measuring the frequencies is shown 
in Fig. 2. A pulsed supersonic beam of CH radicals, with a repetition rate of 
10 Hz, is produced by photo-dissociation of CHBr 3 (96% purity stabilized in 
ethanol) in a 4-bar carrier gas of either He, Ne, Ar or Kr (refs 46,47). The 
measured mean speeds of CH for each of these carrier gases are v 0 = 1710, 810, 
570 and 420ms _1 , and the translational temperature is typically 0.4-0.5 K. 
About 95% of the molecules are formed in the /= 1/2 state. The mixture, formed 
by bubbling the carrier gas through liquid bromoform, expands through the 
1-mm-diameter nozzle of a pulsed solenoid valve into the vacuum chamber. 
The light from the photo-dissociation laser (220 mj, 20 ns, 248 nm) is 
focused in front of the nozzle to a spot 1 mm high (along y) and 4 mm wide 
(along z). 

At z = 86 mm, the molecules pass through a 2-mm-diameter skimmer and 
enter the second chamber, where the pressure is below 10 _ 7 mbar. They pass 
through the state selector and transmission line resonator described below, 
and are then detected at z = 780 mm by driving the A 2 A(v = 0) <-X 2 Il(v = 0) 
transition with a laser and imaging the resulting fluorescence onto a 
photomultiplier tube. In this way, the time-of-flight profile of each molecular 
pulse is recorded with a temporal resolution of about 5 |is. The probe laser 
beam comes from a frequency- doubled continuous- wave titanium- sapphire laser. 
The beam propagates along x, is linearly polarized along z, has a wavelength 
near 430.15 nm and has a power of 5mW in a rectangular- shaped beam, 4 mm 
high and 1.4 mm wide. We measure the population in the (1/2 _ , F) levels for 
the /= 1/2 measurement and in the (3/2 + , F) levels for the /= 3/2 measurement, 
by driving the transitions to the / = 3/2 and / = 5/2 levels of A 2 A(v = 0, N=2). 
The frequencies of these transitions, which are labelled £ 2 2ff(l/2) and 
_Rnff(3/2), are given in ref. 48. The CH density at the detector is approximately 
10 6 cm~ 3 . 

The state selector is situated at z = 241 mm. For the /= 1/2 measurement it 
is a laser beam of the same frequency, polarization and propagation direction 
as the detection laser. This beam depletes the population of the negative parity 
component. For the J = 3/2 measurement, we drive population from (1/2 ~, F) to 
the (3/2 + , F') state, using approximately 10 uW of radiation near 533 GHz. This 
mm-wave radiation is the amplified 54th harmonic of a frequency synthesizer. The 
radiation is collimated by a teflon lens of 30 mm focal length, giving a 10-mm- 
diameter beam, which propagates along x and is polarized along y. The efficiency of 
the population transfer is about 50%. 

The microwave transmission line resonator is built from a pair of 30-mm-wide 
parallel copper plates separated by 5 mm along y. At the A-doublet frequencies of 
interest, the transmission line can only support a TEM 00 mode polarized along y. 
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Microwave radiation is launched into one end from a semi-rigid coaxial cable. 
The other end is not terminated and so a resonator with a quality factor of 
approximately 10 is formed. For each frequency component measured, the length 
was cut so that the transmission line had a resonance at that frequency. The 
length was approximately 480 mm for the /= 1/2 measurement and 440 mm for 
the /=3/2 measurement. The microwave radiation is generated by a synthesizer 
phase locked to a GPS frequency reference whose fractional accuracy is better 
than 10 ~ 13 . A fast, high isolation switch controls the microwave pulse duration. 
The transmission line is housed inside a mu-metal magnetic shield. For the 
7=1/2 measurement a single shield is used, whereas for the 7=3/2 measurement 
we use a pair of nested shields. The outer shield extended into the region 
of the state selector, and a pair of holes cut into this shield allowed the mm-wave 
beam to pass through. Coils inside the shields allow us to apply well-controlled 
magnetic fields. 

For each transition, we first make an approximate measurement of the 
transition frequency by applying a single microwave pulse to the molecules and 
measuring the population as a function of the frequency. Owing to the two 
counter-propagating waves, we observe two peaks separated by twice the Doppler 
shift. Their relative amplitudes indicate that A ^0.1. The mean frequency of these 
peaks is the transition frequency. With the frequency fixed at this value, we map 
out the standing microwave field by driving the transition with a 15-u.s pulse and 
scanning the time when it is applied. The times at which the resonance is largest are 
the times when the molecules pass each of the antinodes. Next, we apply a short 
pulse when the molecules are at an antinode, and we scan the microwave power. 
We observe Rabi oscillations in the population, and thus determine the power 
needed to drive a TT-pulse. We do this for each of the antinodes. With the 
parameters now fixed, we measure the transition frequency using the Ramsey 
method of separated oscillating fields 39 , as described in the main text. Compared 
with the continuous wave excitation, this method has a higher resolution and is 
relatively immune to the broadening and distortion of the lineshape caused by non- 
uniform fields. 
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